PARTIAL TRANSPOSE OF RANDOM QUANTUM STATES, 
WEINGARTEN CALCULUS AND MEANDERS 



MOTOHISA FUKUDA AND PIOTR SNIADY 

Abstract. We investigate the asymptotic behavior of the empirical 
eigenvalues distribution of the partial transpose of a random quantum 
state. The limiting distribution was previously investigated via Wishart 
random matrices and shown to be the semicircular distribution or the 
free difference of two free Poisson distributions, depending on how di- 
mensions of the concerned spaces grow. We use instead Weingarten cal- 
culus to study this problem and find three natural regimes in terms of 
geodesies on the permutation groups. Two of them correspond to the 
above two cases; the third one turns out to be a new matrix model for the 
meander polynomials. Moreover, we prove the convergence to the semi- 
circular distribution together with its extreme eigenvalues under weaker 
assumptions, and show large deviation bound for the latter. 



1. Introduction 

In this paper, we investigate asymptotic behavior of the empirical eigen- 
values distribution of the partial transpose of the random quantum state 
(positive semidefinite Hermitian matrix of trace one) generated by the Haar 
measure on the unitary group. This problem originated from the field of 
quantum information theory in relation to detecting entanglement in a bi- 
partite system. Non-entangled states, called separable states, are necessarily 
positive semidefinite after partial transpose [Per96|, where the latter prop- 
erty is called Positive Partial Transpose, abbreviated as PPT. The converse 
statement is not true except for bipartite states on C 2 <8> C 2 and C 2 <g> C 3 
[HHH96J. Here, the partial transpose is made by writing the bipartite matrix 
as a Kronecker product and transposing each block. Therefore, the generic 
eigenvalue distribution of the partial transpose of a random quantum state 
is interesting and especially the behavior of the minimal eigenvalue is im- 
portant. 

Mathematically, we investigate the following problem. Take three com- 
plex vector spaces C L , C m and C n with I, m, n e N and define p = Tr c i \v) (v 
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for a uniformly random unit vector \v) in the product space C' (g> C m (g> C n . 
Space C z is called the environment, space C m ® C" is called the system, 
individual spaces C m and C n correspond to individual parts of the bipar- 



tite system. We will present the details of this construction in Section 2.4 



This way of in ducing measure on mixed quantum states was investigated in 
[ZS01 , BZ06J. Then, we study the asymptotic behavior of the eigenvalues 
of its partial transpose p r where the transpose acts only on the space C n . 

Although quantum states correspond to Hermitian matrices of trace one, 
Aubrun in HAubl2H used the normal Wishart matrix model to represent a 
random quantum state. The trace of such a Wishart matrix is a random vari- 
able which converges to one only asymptotically. Aubrun showed that the 
empirical eigenvalues distribution converges to the semicircular distribution 
as the dimension of the spaces grow in such a way that / oc mn and man. 
Later, Banica and Nechita [B N12al[BN12bll showed that the limiting distri- 
bution is the free difference of two free Poisson distributions via the Wishart 
matrix model in the regime where the dimension m of one of the parts of 
the system is fixed and the dimension / of environment and the dimension 
n of the other system grow proportionally. 

By contrast, we look into this problem more directly by using Weingarten 
calculus to avoid Wishart matrices and make direct calculations on random 
quantum states. Then, we calculate the moments to find three natural types 
of geodesies which yield interesting distributions. Two of them correspond 
to the above mentioned cases ( HAubl2H and HBN12a[ lBN12b!0 . and the re- 
maining case turns out to be related to the meander polynomials [DFG G971I . 

Our paper is organized as follows. After explaining necessary mathemat- 
ical techniques (Weingarten calculus and free probability) and our precise 
mathematical model in Section [2| we analyze the regime, where I, m, n 
grow such that I oc mn in Section [3} In HAubl2l one requires m oc n but 



we drop this condition to get the same limiting measure in Section 3.1 al 



though we need some weak conditions to show the convergence of extreme 



eigenvalues and their large deviation in Sections 3.2 and 3.4 Then, it is 



shown in Section [4] that our random matrix model yields the meander poly- 
nomials. The connection to free Poisson distribution is presented in Section 
[5] Section[6]contains the concluding remarks. 



2. Preliminaries 

2.1. Free probability, noncrossing partitions and permutations. In this 
paper we will use some basic notions from free probability theory. A good 
treatment of this topic is given in the book [NS06J. We will recall briefly 
the most important notions. 
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2.1.1, Noncrossing partitions. A noncrossing partition r is a partition of 
the set [n] := {1, . . . , n} with a property that if a < b < c < d G [n] are 
such that a, c belong to the same block of r and b, d belong to the same block 
of r then a, b, c, d belong all to the same block of r. Noncrossing partitions 
can be represented graphically as noncrossing connections between points 
arranged on a circle, see Figure [T] 

The set of noncrossing partitions of [n] will be denoted by NC(n). We 
also use the notation 

NC il ,.„ )il (n) := {rGNC(n) : |c| G Vc G r} 

for the set of noncrossing partitions of [n] with a restriction on sizes of the 
blocks. 

2.1.2. Permutations. We denote the permutation group of n elements by 
S n . For a permutation a G S n we define #a to be the number of cycles 
in a and define its length \a\ as the minimum number of factors necessary 
to write a as a product of transpositions. Then the distance in the Cay ley 
graph 

dist(a, (3) := 

is a metric on S n and 

#ct = n — \a\ 

holds for all a G S n . 

For given permutations a, 7 G S„ we define the geodesic a — > 7 as the 
set of all permutations which are on the geodesies between a and 7: 

a — > 7 := |/3 G £„ : dist(a, /3) + dist(/3, 7) = dist(a, 7)}. 

If /3 belongs to this geodesic, we denote it by a — > (3 — > 7. 
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Figure 2. (a) Graphical representation of the permutation 
t(r) = (1) (236) (45) from the geodesic id ->• tt = (12 . . . 6) 
corresponding to the noncrossing partition r from Figure [T] 
(b) Graphical representation of the permutation (t(r)) 1 = 
(T) (632) (54) from the geodesic id -> vr" 1 = (654... 1) 
corresponding to the noncrossing partition r from Figure [TJ 



2. 1 .3. Noncrossing partition and permutations. We will recall now the con- 
struction of Biane HBia97l . We consider the canonical full cycle 



7T 



1,2, 



,n) E S n . 



For a given partition r of [n] and i E [n] we define (t(r))(i) E [n] to be 
the element in the same block of r as i which is after i (with respect to the 

2a Formally, (t(r)) (i) is the first element of the 
Jjt(i)) , • • • which belongs to the same block of r 



cyclic order), see Figure 
sequence 7r(z 



7T 



as i. It is easy to check that t(r) E S ri 

It is easy to check that (i(r)) X (i) E [n] is the element in the same block 



of t as i which is before i (with respect to the cyclic order), see Figure 2b 



The relationship between noncrossing partitions and geodesies in the 
Cayley graph of the symmetric groups is given by the following result. 
Comparison of Figure [T] and Figure [2] is probably the best way to give an 
intuitive meaning to this relationship. 

Lemma 2.1. 

a) The map r h-» t{r) is a bijective correspondence between NC(n) 
and the geodesic id — > tt. 

b) The map t4 (^( r )) is a bijective correspondence between ~NC(n) 
and the geodesic id — > n^ 1 . 
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c) The map r \-t t{r) = (t(r)) 1 is a bijective correspondence be- 
tween NCi ) 2(n) and the intersection of geodesies (id — > n) fl 
(id-^vr- 1 ). 

Proof. Part |g)| was proved by Biane [Bia97j. Part K)] follows in a similar 



way be replacing permutation 7r by tt 

NCi i 2(n) then each cj 

it follows that 



We will prove now part c) If r G NCi^n) then each cycle of t(r) has 

-l 



length 1 or 2, thus t{r) = (i(r)) . From part a) and b) 
t(r) = (t(r)) 1 belongs to each of the two geodesies. 

In order to show surjectivity, let a E ( id — > 7r) fl ( id — > 7r _1 ). From 

part a) we know that there exists r G NC(n) such that t(r) = a. From part 

-l 



b) we know that there exists r' E NC(n) such that (t(V)) = a. Since 
permutations t(r) and t(r') differ just by orientation of the cycles, it follows 
that r = t' . Thus t{r) = a = (t(r)) 1 ; it follows that a is an involution, 
therefore each block of r consists of 1 or 2 elements. Therefore we showed 
existence of r G NCi )2 (n) such that t(r) = a. □ 

2. 1.4. Genus functions. It will be convenient to consider the following two 
non-negative, integer functions on S n given by: 

2g^\a) := dist(id, a) + dist(a, 7r _1 ) — dist(id, 7r _1 ), 
2g^\a) := dist(id, a) + dist(o;, 7r) — dist(id, 7r), 

for a E S n . They are called genus functions; they measure how the paths 
via a are longer than the geodesic distance between id and n~ l or n. 

2. 1 .5. Free cumulants. For a probability measure /i on the real line R hav- 
ing all moments finite we consider its/ree cumulants (^ p (/i)) p l 2 given 
by the following implicit relationship with the moments: 



(2.1) m p (n) := J x p dfi(x) = ^ k\ b \ (/i) . 

For example, 



TSNC(p) b£r 



mi —ki, 

m 2 =k 2 + k\, 

TTi3 =kz + 3kxk 2 + k\, 

m A =k 4 + Ak x k 3 + 2k\ + <ok\k 2 + k\. 

Free cumulants are a fundamental tool of the combinatorial approach to free 
probability theory HNS06L 
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2.1.6. Semicircular distribution. The semicircular distribution with mean 
M and standard variation a, which will be denoted by SCm,<t, has the fol- 
lowing density: 

dSCM >° = — J^a 2 -(x- M) 2 for \x - Ml < 2a. 
ax 2ira 2 

This measure has a compact support [M — 2a,M + 2a] . The free cumulants 
of this measure are given by 

?M ifp=l, 
a 2 ifp = 2, 
otherwise. 



In the special case when M — 1, the moment-cumulant formula ( |2.1[ ) im- 
plies that the moments of this measure are given by: 

m p (SC liCT )= Yl a 2 ^ ofblocksonen ^ 2iaT \ 

reNCi )2 (p) 

Note that any r G NCi j2 (p) can be identified with an involution r := t(r) = 
(t(r)) 1 e S p . By this identification, we have 

2 | T | 

tSNCi, 2 (p) 



(2.2) m p (SC 1>CT )= ^ o- 



2.1.7. Free Poisson distribution. Let A > and ael, Then, the free Pois- 
son distribution with rate A and jump-size a is defined to have the following 
probability density v on K: 



(2.3) 2A Q 



;i-A)5 + Az> if < A < 1, 
v if 1 < A < oo. 



Here, u x , a is the measure supported on the interval [ck(1— Va) 2 ,a(l + %/A) 2 ] 
with density: 

*a,«(*) = ^^^-(t-^a + A)) 2 ^. 

Importantly, free cumulants of this distribution are particularly simple: k n = 

Xa n . 

When a = 1, the free Poisson distribution is called in particular Marcenko- 
Pastur distribution (with variance 1 and parameter 1/A). 
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Figure 3. Meander of order n = 4 with 
k = 2 connected components corresponding to 

a x = {{1,2}, {3, 4}, {5, 8}, {6, 7}} and a 2 = 
{{1,6},{2,5},{3,4},{7,8}}. 

2.2. Meanders. A meander of order n is defined as follows. Suppose we 
have an infinite straight river and 2n bridges over the river. Then, a mean- 
der is a collection of closed self-avoiding and noncrossing connected roads 
passing through all of the bridges; in other words, a meander of order n con- 
sists of a number of loops crossing a straight line at 2n points, see Figure [3} 

(k) 

For each k — 1, 2, . . . n, we define M„ to be the number of nonequivalent 
meanders with k connected components. Another formulation of this object 
is via noncrossing pair-partitions. If one crosses a bridge then next this per- 
son must cross another bridge to come back to the original side of the river; 
these choices can be represented by an element of NC 2 (2n) for each side of 
the river. Therefore, equivalent classes of meanders are represented by the 
elements of NC 2 (2n) x NC 2 (2n), i.e. pairs of noncrossing pair-partitions. 
Then, the meander polynomial M n (x) is defined to be 



M n (x) := J2 xkM n k) - 



k=l 



This polynomial has been extensively investigated and some matrix models 
were found in [DFGG97]. However, we think that our model can be the 



simplest matrix model for the meander polynomials (see Theorem 4.2 ). 



2.3. Weingarten calculus. A graphical version of Weingarten calculus in- 
troduced in HCN10H enables us to transform the algebraic calculus into a 
graphical one. This method turns out to be very useful, see for example 
HCFN12H . We briefly review the graphical method in this section. 

We have the following formula for integrals of polynomial functions in 
entries over the unitary group U(n) with respect to the Haar measure: 

Theorem 2.2 (Weingarten [WeiTU). Let n be a positive integer and i = 
(ii, . . . , i p ), i' = . . ■ , i'p), j = (ji, • • • ,j P ), f = (ii, ■ ■ ■ J'p) be p-tuples 
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of elements of [n]. Then 

(2-4) / U iljl -.-U ipjp U^...U^,dU 

Ju(n) 



Ifp ^ p' then 



E5i,i' ■ ■ ■ 5i i> 5i,j' ■ - • Si ,•' Wgfn, da 1 " 



u ixh ■ ■ ■ u i P j P Ui'Ji ■ ■ ■ u t,?„ du = °- 



p' p' 



>U{n) 

In the above th eorem , Wg(-, ■) is called Weingarten function and it is 



deeply studied in IICS06I . However, in the paper we only use the following 
property: 

Lemma 2.3 (Collins [Col03]). For any value ofn 

P-i . 

]TWg(n,a) = n— • 

a&Sp i=0 

We overview the graphical calculus; we refer the reader to HCN10B for 
detailed explanations. The graphical calculus consists of boxes and wires. 
A box represents a vector, a dual vector or a more general tensor (such 
as a linear map), which is essentially a linear combination of products of 
vectors and dual vectors. A wire is a generalization multiplication, trace, 
and similar operations. For example, we view the following two operations 
as identical in the sense that both are described by a kind of delta function: 

multiplication : AB = '^^a i jbki \i)(j\k)(l\, 



ijkl 



trace : Ti A = VVjli) 



a>ij- 



Here, A = a ij I*) 01 and B = Y,ki hi\k}(l\. In other words, by using 
wires, we can indicate how the structure of the contractions of the tensors 
(or, viewed alternatively, indices). 

In Figure |4} white circles and black circles represent primal and dual parts 
of linear maps, respectively. If a linear map acts on a product space, we can 
use several symbols to represent different spaces. Wires themselves can be 
interpreted as a generalization of the identity map or the Bell state, etc., but 
we do not go further in this paper. 



The right hand side of ( |2.4[ ) involves only delta-functions and the Wein- 
garten function. This means that through an integration over the unitary 
group we can eliminate the boxes of U and U and reconnect wires which are 
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B 




(A) 



(B) 



Figure 4. Comparison of (a) multiplication and (b) 



trace. 



attached to those boxes according to the delta-functions on the right hand 



side of ( |2.4| ). This process for fixed permutations a, (3 is called removal and 
the whole process (which sums over all choices of a, /3 the corresponding 
new graphs) is called graph expansion. 

2.4. Our model. Suppose we have three complex vector spaces C l , C m 
and C n with l,m,n G N and take the uniformly random unit vector \v) in 
the product space C' <8> C m ® C n . Then, the corresponding random pure state 
\v)(v\on C'®C m <g)C n induces a random mixed state p = VV* on C m ® C n 
by tracing out the space C'. Here, we use the usual isomorphism |t>) h->- V 
between (g) C m <8> C n and Ai mn ,i (C) . We study the asymptotic behavior of 
the eigenvalues of its partial transpose p r with the transpose acting only on 
the space C". More precisely, let {At}™} be the eigenvalues of the rescaled 
random matrix mnp r ; we define the corresponding empirical eigenvalues 
distribution 

1 



mn 



our goal is to find the limiting measure in the sense of weak convergence. 
Note that this scaling is used in Section[3]and some other scalings are chosen 
in the following sections. 

First, we transform p into the graphical representation. The uniform ran- 
dom unit vector in C^®C m ®C n can be realized concretely as U \ei)\e m )\e n ) 
where U is chosen randomly from the unitary group U(lmn) with respect 
to the Haar measure and |ej), \e m ) and \e n ) are fixed unit vectors in the 
respective spaces. Hence, algebraically, 

p = Tr c * \U(E l ®E m ®E n )U*} 

is a random mixed state on C m <8> C™, which we are interested in. Here, 
E\ = \e\)(e\\ for A = l,m,m are appropriate rank one pro jections. Thi s 
model of random mixed quantum states was investigated in [|ZS0ll lBZ06l. 
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Figure 5 . A graphical representation of p. 

The graphical representation of p is shown in Figure [5j where circles 
correspond to C l , squares to C m and diamonds to C n . For simplicity, we 
omit those symbols for E\ matrices (A = /, m, n). Importantly, we use U 
instead of U* since the Weingarten calculus is written in terms of U and U. 
As a consequence, we must switch prime and dual parts of the matrix; this 
explains why the white symbols are outside and the black ones are inside. 

2.5. Moments of p r . For an integer p > 1 we consider the random variable 

Z ( n p) := — Tr [(mnp r ) P ] = [ x p dp mnp v 
mn J 

which is just the appropriate moment of the empirical eigenvalues distribu- 
tion p. mrip r. The following theorem will be our main tool. 

Theorem 2.4. For arbitrary integer p > 1 the expected value of the corre- 
sponding moment ofmnp r is given by 

(2.5) EZ^ p) = F{lmn,p) Z ~ H ™ P-1 " M n p - l -^ la \ 

,~ -w, x fmn\\ a \ 9 „(i) 9 (2) 

(2.6) =J(l mn ,p)^ - m' 29 " n~ 29p , 



V I ) 

aeSp 



where := gp\a) were defined in Section 



2.1.4 



and 



(2.7) F{D,p) := D p Wg( J D, 7 ) = JJ — = 1 + (-D -1 ) . 

-yeSp i=0 

Proof. Before calculating ETr(p r ) p we must understand (p r ) p , a part of 
which is drawn in Figure |6] Notice that the transpose on the space C n im- 
plies that the connections between white diamonds have a different structure 
than the connections between white squares. 
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u 



Figure 6. A part of a graphical representation of (p r ) p . 



Fixed permutations a, (3 E S p (as in Theorem 2.2) reconnect white and 
black symbols respectively via the removal. As a result, the white symbols 
give the following factors: 

• white circles: l#( a \ 

• white squares: m*^ a \ 

• white diamonds: la \ 

On the other hand, black symbols always give contribution equal to 1. 
Indeed, since E x (for A = l,m,n) is a rank-one projection, any resulting 
loops with boxes, called necklaces in HCFN121 . always give contribution 
equal to 1. 

Hence, 

E Tr[(/) P ] = V / #(Q) m #( ™ ) n #( ^ la) Wg(/mn,a- 1 /3) 

CeW(imn) ^ — 4 



F(lmn,p) r '° 

a£S p 



I m -\ na \ n -\ n 1q: I. 



and the result follows by simple algebraic manipulations. □ 



Interestingly, the above formula ( |2.5[ ) gives three regimes which give in 



teresting limiting measures. The first case, which we investigate in Section 
[3} is that / oc mn, where three permutations a, rra, n~ l a interact. The fol- 
lowing sections treat the other two cases when / or m are fixed so that just 
two of the permutations a, na, n^ 1 a interact. 

3. The case when dimensions of environment and both 
system parts are large 

In this section we investigate the regime where l,m,n — > 00. Aubrun 
HAubl2l investigated the case where / oc d 2 and m, n oc d as d — > 00 via 
Wishart random matrices model and showed the results which correspond 
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to Theorem 3.1 and Theorem 3.2 However, we instead investigate asymp- 
totic behavior of the eigenvalues of p r in a more general setting by Wein- 



garten calculus. We will show the following results. Theorem 3.1 shows 
that mnp r has the limiting distribution as long as the dimension mn of the 
quantum system C m ® C " is proportional to the dimension / of the environ- 
ment C'. Then, Theorem 3.2 studies the behavior of the smallest and largest 
eigenvalues. They turn out to converge to the two corresponding edges of 
the support of the limiting density unless m and n grow too differently. Fi- 
nally, large deviation property of the extreme eigenvalues is investigated in 



Theorem [3T4| under an additional assumption that m and n are not too close 
to each other. 

3.1. Limiting eigenvalues. We analyze the limit of the empirical eigen- 
values distribution of mnp r when / oc mn where l,m,n — > oo. In the 
following, we assume without loss of generality that m > n. 

Theorem 3.1. Suppose — > a with < a < oo and m > n. Then, as 
n — )■ oo, the empirical measure ofmnp T converges weakly almost surely to 
the semicircle distribution SCi^. Here, we think of m = m n and I = l n as 
sequences which implicitly depend on n. 

Before presenting the proof we remark that in the case a = the limit 
distribution SCi o = <$i becomes a delta measure. 



Proof of Theorem 3.1 Non-zero contribution to ( |2.6[ ) in the limit is given 
only by the summands for which a is on the following two geodesies: 

id — t- a — > 7T _1 ; id — > a — > n. 



We apply Lemma 2.1 :) , it follows that 

a e NC 1)2 (p). 

Hence, we have 



(3.1) 
with 



lim EZi p) 



E 



m p (SCi i(r ) 



aGNCi, 2 (p) 



a~ 



where we used ( |2.2[ ). Thus, we proved the convergence in (expected) mo- 
ments. 

To prove almost sure convergence, we will show later that 

oo oo 

(3.2) ][>rZ<f) = ;>> \{Z^f] - (EZ<f)) 2 < oo 



71=1 



n=l 
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for each p E N. This result via standard arguments (involving Markov in- 
fo) 

equality and Borel-Cantelli lemma) would imply that Zn converges almost 
surely to the appropriate moment of the semicircle distribution. The latter 
distribution is uniquely determined by its moments, so convergence of mea- 
sures in the sense of moments implies the convergence in the weak sense; 
this would finish the proof. 

It remains to show that ( |3.2| ) indeed holds true; we shall do it in the 



following. A more careful analysis of ( |3 . 1 [ ) shows that every term which 
converges to zero is in fact at most O (n~ 2 ) thus 

EZ?>= £ (™) w +0 („->). 

aeNCi l2 (p) 

We consider permutation 

7T 2 := (l,2,...,p)(p + l,p + 2,...,2p) E S 2p . 



One can easily check that E [Z^ J is equal to the right-hand side of ( |2.5[ ) 
and ( |2.6[ ) in which the summation over S p is replaced by summation over 
S 2p and permutation n is replaced by n 2 . In an analogous way it follows 
that 

*(*?>)' =( £ (™) w ) >+ °(»-) 

Va6NGi, 2 (p) / 



2 



which finishes the proof of (3.2 ). □ 



Theorem 3.1 shows that the limiting empirical distribution has the com- 
pact support on the interval [1 — 2y/a, 1 + 2y/a\. However, this does not 
necessarily mean that the minimum and maximum eigenvalues converge 
to the boundaries of this interval. We analyze the convergence of extreme 
eigenvalues in the next section. 

3.2. Behavior of extreme eigenvalues. In this section we analyze the be- 



havior of minimum and maximum eigenvalues. Theorem 3.2 below shows 



that in the regime of Theorem |37T the minimum and the maximum eigenval- 



ues of mnp r actually converge respectively to 1 — 2y/a and 1 + 2y/a under 
an additional condition on growth of m and n. 

3.2.1. Convergence of extreme eigenvalues. 



Theorem 3.2. Let assumptions of Theorem 3. 1 hold true with an additional 
condition that logm = o(n 1 / 6 ). Then the extreme eigenvalues of mnp T 
converge to 1 ± 2^fa almost surely. 
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Proof. The difficult part of the theorem is to show that 
limsup \\mnp v — 1]]^ < 2\fa 

n 

holds almost surely. We will do it in the following. 

Let {p n }n be a sequence of even numbers such that p n is the largest even 
number such that 2p 12 max {l, ^y 1 } < n 2 ; we write p = p n for simplicity. 
Note that then logm = o(p). 

Hence we have 

E \\mnp r - < E \\mnp r - I\\ p p = 

E Tr [{mnp v - If] < 2mnp 5 (2y/a + o(l)) P . 

Here, the last inequality follows from Lemma[33]below. Therefore, Markov 
inequality implies that for any e > 



p 

< oo. 



Y Pr { \\mnp r - 1^ > 2^~a + e} < V 2mnp 5 ( 2 V E + °( 1 ) ) 

' V 2 V a + e / 

n n x v ' 

Thus Borel-Cantelli lemma finishes the proof. □ 
Lemma 3.3. If m > n and 2p 12 max{l, a} < n 2 , then 

Vap\ P 



ETr [(mnp r - l) P ] < 2mnp 5 2y/a + 



mn 



where a := m f k . 



Before presenting the formal proof let us make an intuitive remark on this 
phenomenon. Since we now know that mnp T obeys the shifted semicircle 
law in the limit, we must more or less have the following: 

— ETr [{mnp T - l) P ] « (^) p Cat p/2 « (2^) p , 



where used the fact that Catalan numbers 

Catfc = kTJ (^) = / ^SC ,i(x) 
are the moments of the centered semicircular distribution. 
Proof of Lemma \373\ First, we expand 

— ETr \(mnp T - l) P ] = V ( —ETr \(mnp r ) k ' 

mn \kj mn L 

fc=o v 7 v 

(*) 

Here, ( |2.6| ) implies that 

(*) = F{mnl, k) y a H m~ 2ff * V 2 ** . 

ae5 fc 
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Next, we decompose Sk in the following way: 

Sk = Sk,o U • - • U Sk,k, 

where 

Sk,t — { a G Sk '■ # (non-fixed points of a) —t}. 

For example, Sk$ = {id} and Sk.i = 0, etc. We also define another nota- 
tion: 

St = St,t c St. 

For a given d G Sj^t we consider its support, i.e. the set of its non-fixed 
points and the restriction of a to the support. This support has t elements; 
by a relabeling of these elements, the restriction can be identified with a 
permutation a E S t . If the relabeling is chosen to be order-preserving, one 
can easily check (for example, by removing the fixed points one by one) 
that 

|&| = |a|, 
# (an k ) = # (ant) , 
# (cT 1 ^) = # [pc-^t) , 
thus the genus functions on a and a are the same: 



fori = 1,2. 
Hence, 

— E [(mn/ - 
mil 



*) n -2si 2) («) 



a |a| m -2 9t (1) ( a ) n -2« y f ) (a) 



ST (P\{-iy- k F{mnl,k) a H m- 2 ^^ a 

fc=0 ^ ' a£S k 

t(?)(-ir^(-<*)E(!)E 

fc=0 V 7 t=0 V 7 a(zSt 

Y E oHm"^W n "^W £ (^) ( k ) (-ir k F(mnl, k) 
t=0 ae s t fc=t V / V / 



(*) 



Here, the last inequality comes from Lemma A. 1 
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Next, define a function h t : 
(3.3) ht = gl 1] + gl 2) . 



For this function we claim that for a G S, 



t . . t , . . 

Here, the first inequality comes from the definition of S t and the second one 
is proved similarly as in ( |A.4| ). Then, we have 

(*)<E« M »-* W < E E f" 8 ^ 1 )" 

a£St h>0aGT tM V 7 

E, t 19h , r t / maxjl, a} 

h>0 



\ n 2 
p 12 max{l, a} 



n 2 



h>0 

< 2p B (2y/a) t . 

Here, the first inequality comes from m > n, the third one is implied by 
Lemma [A.21 

Thus we showed that 

^E[(mnp r -/) P ] < Vf^ f^Y~V (2^. 
ran \t J \ van J 

Finally, the binomial formula gives our desired result. □ 

In the following section, we specify how fast those extreme eigenvalues 



can converge as is in Theorem 3.4 



3.2.2. Speed of convergence. Below, we show a large deviation bound for 
the extreme eigenvalues by using an extended version of Levy's lemma in- 
troduced in [ASW11]. Again, in order to get a nice bound we need further 
additional condition on the proportion — v a. We prove the result only 
for the minimum eigenvalue since the proof for the maximum eigenvalue is 
almost the same. 



3.2 



B.2 



assume in addition that 
Then, we have the fallow- 



Theorem 3.4. In the regime of Theorem 
max{a, 1/a} > C 2 where C is as in Lemma 
ing large deviation bound for the convergence of the minimum eigenvalue 
of mnp v to 1 — 2\fa: there exist three constants < ei, C\ < 1 and l\ € N 
such that 

P r {A m in < 1 — 2yfa — e} < exp{— Cxe 2 /} 
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for all < e < e\ and I > l\. 

Proof. We denote the unit sphere in C d by Sj. Then define a function A m j n : 
<S /mn -> R by: 

Amin(^) = (the minimum eigenvalue of mn(VV*) T ) . 

Here, again we use the identification \v) = V between the two spaces <g> 
C m <g) C n ~ -M mrii /(C). For our purposes, this function must be modified 
similarly as in HASWllH . Let K C S; mn be the compact subset of the whole 
domain of A m i n defined by 

K= Iv e S lmn 

First, we claim that the probability of its complement is exponentially 
small in I: 

(3.4) Pi{VeK c ] < exp ( -cmin{l,a} /). 

nd, we calc 
When a > C 2 we 




B.2 



Here, the constant c is as in Lemma B.2 To verify t his b ound, we calculate 
the quantity + which is involved in Lemma 
set di = / and d 2 = mn so that 

i c i r i i 



V~l V™n V~l VI 
On the other hand, when - > C 2 we let rfi = mn and rf 2 = I and have 

'a — 1 ^ 

1 C7 1 C 1 1 



\A7l Vol Vol 

Next, assume without loss of generality that A min (V) > \ m - m (U) and 
denote the normalized eigenvector for the latter by \u). Then, writing p v = 
VV*, we have 

Amin(^) - A min (£f) < (w|mnpy|w) - (u\mnp V jj\u) 



< 


mn 


\p T v 


-P T u\ 


1 oo 


< 


mn 


Vv 


-P V u\ 


2 




mn\\ 


Pv 


- Pu\\ 


2 


< 


mn( 


\\v\ 


oo T 


u 



)\\V-U\\ 2 

< 4max{a, Va}Vl\\V - U\\ 2 . 

Hence, the Lipschitz constant of A mm on K is upper-bounded by 4 max{a, V®} V~l- 
Finally we extend this restricted function to the whole domain by: 



Amin 



(V) = M K {\ min (W) +4max{a, Va}Vl \\V - W\\ 2 \ for V G K°. 
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This is a modified function of A m i n and different f rom the original only on 
the small domain K c . This implies, via Theorem 3.2 that for each e > 
there exists some In EN such that 

e 



(l - 2y/a) - Median 



A, 



< 



for all I > In. 



On the other hand, by applying Lemma BTT] (Levy's lemma) to Si mn we 
have 



(3.5) Pr | Median 



A, 



r f , A ^n f c (2a/ 2 - 2)6 2 



exp 



c (l - (al)- l )e 
32 max{a, 1} 



To finish the proof, set for example 



'32 max{a, 1} x cmin{l, a} /32ca 



c V c o 

and choose appropriate constants c\ and l\ so that applying the union bound 
method to ( |3.4| ) and ( |3.5[ ) gives the desired result. □ 



3.2.3. Implications for quantum, information theory. As already shown by 
Aubrun HAubl2L a = 1/4 is the critical value. The smaller edge of the 
support of SC11/2 is 1 — 2y/a, which is strictly negative when a > |. In 
this case Theorem 3.1 implies that the limsup of the minimum eigenvalue 
of the matrix p T converges almost surely to a strictly negative number. So, 
our random quantum states are generically not PPT when a > \ . On the 
other hand, when a < |, 1 — is by contrast strictly positive. By The- 
orem |3]2} we know that the probability of quantum states not being PPT is 
exponentially small in n 2 . In other words, our random quantum states are 
typically PPT when a < \ . 



4. The case when the environment space is fixed and its 
connection to meanders 

In this section, we investigate our model when the dimension / of the 
environment space is fixed. Unlike the previous regime, we do not have 
double-geodesies any more. However, interestingly, this regime gives a 
simple matrix models for the meander polynomials. 

4.1. Our model. In this section we investigate the case where 

m 

I = l , > c, 

n 
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for fixed / G N, c > 0, in the limit as n — > oo. As usual, m = m n 
depends implicitly on n. Then, we are interested in the following empirical 
distribution of lmp r : 

(4-1) fJLhnff := 52 

where Aj are the eigenvalues of Zmp r . Here we use a different scaling from 
the one in Section [3] The moments of ( |4.1[ ) can be written as 

— E Tr [(Zmp r ) P l = 

mn ueu(imn) L v y J 

(l + O (^~ 2 )) /pH a l c P-M™l n P-2-|™l-K~ la l_ 



Here, we used (2.5 ) and as before, 



(the power of n) = p — 2 — (\na\ + \n 1 a\) 

■1 if pis odd, 
if p is even. 



< p-2 - |?r 



2| 



Here, as before it — (1, 2, . . . ,p) is the canonical full cycle. This means in 
particular that all the odd moments vanish. When p is even, the bound is 
satiated if and only if a is on the following geodesic: 



7T 1 — > a — > 7T. 



This implies that, for even peN, 
(4.2) lim — E [(Zm/fl = V l* a 



c #(7ra)-l^ 



n-»oo mn UeU(lcn 2 ) 



Unfortunately, the above geodesic is very difficult to understand fully. 
So, we investigate two restricted cases where I = 1 or c = 1. The former 
treats random pure states on the bipartite system, and, interestingly, the 
latter show a connection to the meander polynomials. 

4.2. Random pure states. We start with the case / = 1 of trivial environ- 
ment space C'; this corresponds to p being a random pure state. 

Theorem 4.1. Let I = 1. Then the empirical eigenvalues distribution ( 4.1 ) 
converges almost surely as n oo to 

(4-3) (l - I) S + - c P BVXIX2 

in the weak topology of probability measures. Here, Psy/ XiX 2 zs tne distri- 
bution of By / X 1 X 2 , where the distribution of random variables Xi and X 2 



20 



MOTOHISA FUKUDA AND PIOTR SNIADY 



is the free Poisson distribution v\ fi as in (23), and B takes the value 1 or 



— 1 with probability 1/2, and they are all (classically) independent. 

Before showing the proof let us make two remarks. First, when c < 1 
the coefficient of So from the first summand is negative, but this negativity 
is canceled by the mass from the second summand. Secondly, for I = 1 
operator p r has at most min(m, n) 2 non-zero eigenvalues (this follows from 
analysis of Schmidt coefficients of ip, where p = see HAubl2l proof 



of Proposition 9.1]) which explains the atom at zero at the distribution (4.3 ) 



Proof of Theorem 4.1 For a technical reason, we consider the following 



rescaled empirical distribution: 

mn 

(4.4) — V A, 

mn 



Let p = 2q; we consider permutation a £ S p which contributes to ( |4.2[ ) so 
that permutation r := na belongs to the following geodesic: 

id-+r->7r 2 = (l,3,...,p-l)(2,4,...,p). 

Then the limit moment of ( |4.4[ ) is equal to c x ( |4.2[ ), which can be calculated 
as follows: 



c x (Q = 




2 



id— >T— >TT 2 

Moreover, 

E < #, = E IL 

f eNC(g) f eNC(<j) ver 

coincides with the appropriate moment of the free Poisson distribution v c> \ 
(with rate c and jump size 1). Hence, for the even moments, square root of 
the product of two classically independent free Poisson random variables 
gives the right moments. However, since all the odd moments vanish, we 
must add the factor B to recover our desired moments to get the limit dis- 



tribution of ( |4.4[ ). After rescaling back this distribution, the additional atom 
(1 — -)8 does not change the moments m p of the measure for p = 1,2,... 
but takes care of the correct value of the moment m (the total mass of the 
measure). This shows convergence in expected moments. 

Almost sure convergence can be proven similarly as in the proof of The- 
orem 3. 1 Since the limit is compactly supported, the converges in moments 
implies convergence in the weak sense. □ 

In HAubl2L an unpublished result when I = c = 1 was quoted and it 
reads that the limiting distribution is the product of two independent semi- 
circular distribution SC 1 (with mean and variance 1). Since if Y gives 
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(A) 



(B) 



Figure 7. (a) graphical represent ation of a noncrossing 
partition {{1,3, 4}{2}} G NC(4), |(b)| the corresponding 



noncrossing pair-partition from NC2(8). 

a semicircle law SCo,i then Y 2 has the Poisson distribution v\,\, this is a 
special case of ours. 

4.3. Meander polynomials with our model. Next, we consider the case 
where c = 1, where our model gives the meander polynomials. 



Theorem 4.2. If — — > c = 1 and q G N ?/zen Z7ze 2g-?/z moment of /i; mp r 
converges as n — > oo to ?/ze meander polynomial M q {l). 



Proof. First, since c 
(4.5) 



(ri0T 2 )] 



n,T2 



Here, Ti G NC({1, 3, . . . ,p - 1}) and r 2 G NC({2, 4, . . . ,p}). 

Next, we recall the well-known bijection between NC(g) and NC2(2g). 
We represent a noncrossing partition as in Figure [Taj We add two points i_ 
and i + for both sides of each % G [q], left and right respectively. Then we 
connect i + and j_ if a(i) = j. The example is drawn in Figure 7b where 



we also use arrows to show the action of the permutation a. This procedure 
is commonly known as fattening. 

Finally, we will calculate #[7t -1 (ti © r 2 )]. To understand the loop- 
structure of vr _1 (ri © r 2 ) we note that n and r 2 act in turn; n acts on odd 
numbers and r 2 even numbers, but 7r _1 switches parities. For this reason, 
we suggest the following graphical representation, see Figure [8] 

(1) We draw two parallel horizontal straight lines with odd-numbered 
points on the upper line and even-numbered points on the lower line. 

(2) Draw the graphical representation for t\ above the upper line and 
the graphical representation for r 2 below the lower line. 

(3) Identify respectively (2z + 1)_ and (2z) + , and then (2z)_ and (2i — 

!)+■ 
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Figure 8. Graphical representation of noncrossing parti- 
tions n = (1, 3, 5) (7) and r 2 = (2, 8) (4, 6) 



Here, note that the last step corresponds to the action of 7r _1 . An example 
with is drawn in Figure [8] Since we can think of T\ and r 2 as elements of 
NC 2 (2g), identifying l_x and 2q + reduces ( |4.5p into the meander polyno- 
mials M q (l). □ 



5. The case when one of the system parts is fixed 

In this section, we review the case where one of the spaces C m or C n is 
fixed. Without loss of generality we may assume that C m is fixed. This case 
was investigated by Banica and Nechita in HBN12airBN12bH via approxima- 
tion by the complex Wishart matrix and they proved that the limiting mea- 
sure is the free difference of two free Poisson distributions. We will present 
a new proof of this result. 

Suppose 

1 u 

> b, m = m , 

n 

where b > and m G N are fixed constants. Then, the restatement of 
nBN12allBNT2bll in our setting is: 

Theorem 5.1 (Banica, Nechita). The empirical distribution ofmlp T con- 
verges weakly, as n — > oo, almost surely to the probability measure of the 
free difference of free Poisson distributions with parameters b(m ± l)/2. 
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Proof. First, formula ( |2.5[ ) gives 

— E Tr\(mlp r ) p ] = 

mn U&U(lmn) 



(l + 0(n~ 2 )) J2 bP ( b 



n) H m P-i-M n P-i-k H 



Then, the power of n is bounded as: 

p — 1 — (|a| + |a _1 7r|) < p — 1 — \ir\ = 0. 
This implies the following geodesic formula for the limit: 



lim — E Tr[(m/p r ) p ] = V b*' 

n-»oo mn U&Ailmn) 

id— >a— >tt 



« m #(™)-l 



(5.1) = ^ 6 #c 

aSNC(p) 

Here, for the second equality we used l + e(a) = #(ira), which was proven 
in HBN12a[ Lemma 2.1], where e(-) is the number of blocks whose size is 
even. The latter formula gives the free cumulants of the limiting measure as 



(5.2) k p 
If we set 



bm if p is even, 
b if p is odd. 



b(m + 1) b(m — 1) 



x = — - — ; y 



2 2 
and define X (resp. Y) to be a random variable with free Poisson distribu- 
tion u X; i (resp. fiy t {) then, if X and K are free then the cumulants of X — Y 
are given by 

I o it p is odd. 

so they coincide with the cumulants ( |5.2[ ) of the limiting distribution. In this 
way we showed convergence in expected moments. 

Almost sure convergence can be proven similarly as in the proof of The- 
orem |3.1[ As the limiting distribution is compactly supported, the conver- 



gence in moments implies weak converge. □ 

6. Concluding remark 

In this paper, we investigated asymptotic behavior of eigenvalues of par- 
tial transpose of random quantum states on bipartite systems. We naturally 
picked up three regimes depending on how the concerned spaces grow and 
showed their connections to the semicircular distribution, the free Poisson 
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distribution and the meander polynomials. Other regimes may show other 
interesting behaviors. 



Appendix A. Lemmas for Section [XO] 

Lemma A.l. For arbitrary integers p, t > and D > 2 the following 
bound holds true: 



P -t 



where T is given by ( |2.7[ ). 

Proof. We denote by S the shift operator on functions of a single variable, 
i.e. S : g(-) >->■ g(- + 1) or alternatively (Sg)(x) = g(x + 1). Then, writing 
J-'(k) := J-"(-D, k), the left-hand side of ( |A.1[ ) becomes 

fc>0 v 7 v 7 fc>0 v 7 v 7 



(A.2) = A*> ( 1 ?{■) 

where A = S 1 — 1 denotes the forward difference operator. 

Firstly, we will recover the well known product rule for the finite differ- 
ence A. For arbitrary functions g, h of a single variable we have 

A[gh] = (Sg)(Sh)-gh 

= [(Sg)-g)(Sh)+g[(Sh)-h) 

= (Ag)(Sh)+g(Ah) 

= M(Ag® Sh + g® Ah), 

where M(g ® h) = gh denotes the pointwise multiplication of functions. 
In other words, we showed that 

AM = M(A<g>S + l<g>A). 

It follows that higher powers of the forward difference operator act on prod- 
ucts as follows: 

A p (gh) = (A p M)(g ® h) = ^ (f) x A ^ x S l A p ~ l h. 

i=o ^ 7 

We used the fact that the operators of shift and the forward difference com- 
mute: AS" = SA. 
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By applying this general formula in our particular setup we obtain: 
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(A.2) 



1=0 



P 



=0 



, x A* 
// \t 



X 



(❖) 



(A.3) 



(A p -^)(t), 



where we used the following property of (0): 





(0) = A' 



t-l 



[l = t]. 



Thus our problem is reduced to estimating the quantity ( A p t J z ) (t) appear- 



ing on the right-hand side of ( |A.3[ ). 

If g is an arbitrary function of a single variable then A acts on the product 
g x (S^F) as follows: 

A[g x (S l F)](k) = g(k + l)T{k + i + 1) - g{k)F{k + i) 

= g (k + l)-(l + ±p)g(k)]r(k + i + l) 



A k + i 

A I q 

D 1 " 



D 

(k) x (S m F)(A;) 



where it follows from the definition ([277]) that 

F(k + i) = 
Hence inductively we get 



l + ^jp-)F(k + i + l). 



[A«(g x T)](k) 



A 



k + q- 1 



x (S*F)(fc). 



We are interested in the special case of this formula for g = 1 given by the 
constant function: 



[A^](fc) 



A 



fc + q - 1 



A-*ll 
D 



xf(k + q). 



(*) 
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We will now analyze the structure of expression (4k). We use the short- 
hand notation 

Pi = W) = 

Expression (4) is a product of q factors, each factor being the sum of two 
terms. Let us expand this product; each of 2 q resulting summands is a 
product of the forward difference operators A (let us say that there are r 
factors of this form) and expressions of the form (P?)o<i<q-i (let us say that 
there are q — r factors of this form). Notice that these two expressions do 
not commute so the order of the factors is important. In the following we 
will study in detail expressions of this form. 

Clearly AP; = — thus a straightforward application of the product rule 
shows that , 



I I 

APj(l) • • • Pi(£) = ^2 A P(l) • • • Pj( r ) • • • Pi(£) + 1 



Kr<t 



D ^ 

\<r<i 



P%{\) ' ' ' Pi(r-l)Pi(r+l)+l ' ' ' P(f)+1- 



The right-hand side can be interpreted as follows: we sum over all ways 
of matching the forward difference operator A with one of the factors P( r ) 
on its right; this factor is removed and every term P, on the right of P i(r ) 
should be replaced by Pj+i- This matching has been illustrated graphically 
as an arrow connecting the difference operator A with the factor P i(r ) on 
which it acts. 

This observation can be extended to general products which we consider. 
Namely, we sum over all ways of matching r difference operators A with 
factors (p) in such a way that each operator A is matched with some factor 
Pi which is on its right and each factor P is used at most once. The fac- 
tors (Pi) which are matched should be removed, each unmatched factor Pj 
(there are q — 2r of them) should be replaced by P j+5 where 5 denotes the 
number of factors (Pj) which are matched and are on the left from Pj, thus 
j + S < q — 1. Also, there is additional factor (— j^) ' . 

We can illustrate this by an example for r = 2: one of the summands 
contributing to the product P 2 A P 5 A P 8 P 13 P 17 P 2 \ P25 is given by: 



P 2 A P 5 A P 8 P 13 P17 P21 P 
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I 

+ P 2 A P 5 A P 8 P 13 Pi 7 P 2 i P 25 + • • • = 

• • • + JpPz P5 P% Pis P27 + 
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The above two-step procedure (selecting one of 2 q summands with r fac- 
tors A, then summing over the matching V) can be clearly replaced by sum- 
ming simply over V which should be a partition of the set [q] with r blocks 
of length 2 and q — 2r blocks of length 1 (the places where the difference 
operator A occur are exactly the left elements of two-element blocks of V). 
The number of such partitions V is equal to 



q(q - 1) • • • (? - 2 r + 1) </ 



2r 



2 r r\ 



2 r r\ 



It follows that for k > 
|(*) 

Thus for q = p — t 



< 



E 



1 9 



2r 



0<r<LfJ 



D r 2 r r ! 



k + q 
D 



g-2r 



(A3) 



(A^)(t) 



(*) Hp) 



and 



(A.3) 



< 



E 



2r 



1 p 

W r 2Fr\ 



P \ «-2r 



E 



0<r<LfJ x 7 0<r<L|J 

The sum on the right-hand side is dominated by its last summand multiplied 
by 2, therefore 



(A.3) < 2 



D q l 2 ' 



which finishes the proof. 

Lemma A.2. Let h = h p be a function defined on S p as in (33 ) and 

T Pt h = {« 6 5 P : h(a) = h}. 

Then 



□ 



T p , h \ < ^ 



-lp!2h+5 



Proof. We assume for a while that p is an even number. We define 

k(a) = 2\a\ — p. 

It fulfills the bound 

(A.4) 2h(a) > 2\a\ + |tt 2 | - 2p + 2 = 2\a\ - p = k(a). 

Let cr be an arbitrary permutation and r be a transposition. If \ra\ < \a\, 
we say that a' := to was obtained by a cut. This corresponds to splitting 
one of the cycles of a into two cycles; we say that the cycle was cut. If 
\ra\ > | a\, we say that a' := to was obtained by a gluing. This corresponds 
to merging two of the cycles of a into one cycle. 
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Let a G T Pi h be given. If a has a cycle of length at least 3, we cut it; we 
iterate the procedure until we obtain permutation a' which consists only of 
cycles of length 1 and 2. The number of cuts we need is upper-bounded by 
k(a). To see this, 

(the number of cuts for a) < / ](\c\ — 2) = p — 2#(a) = k(a). 

Here, |c| is the cardinality of a cycle c of the permutation a. Importantly, 

h(a') < h(a). 

Next, we glue fixed points of a' pairwise to get a" which consists of p 
disjoint transpositions; in the language of partition this is a (possibly cross- 
ing) pairing. The number of glues is upper-bounded by k(a)/2; 

(the number of glues) < - ^^(|c| — 2) = -k(a). 



A3 



Each 



2 1 ' 2 

Above, we used the fact that a has no fixed points. 

(2) 

Write g{a) = g p , which is the notation used in Theorem 
operation of gluing can increase the genus at most by 1; furthermore g < h 
so that 

g{a") < g(a') + h{a) < 2h(a). 

To summarize: we constructed a map a i-> a" where a E T P) h, with a 
property that g(a") < 2h and a" is an involution without fixed points. Fur- 
thermore, a can be obtained from a" by multiplying by at most 3h trans- 
positions; this means that the preimage of any a" consists of at most (^j 3h 



elements. Then, by using Lemma A.4 



/ \ 3ft 2ft 
V / g=0 



-1 12/1+3 



which finishes the proof in the case when p is even. 

When p is odd number, we select some cycle consisting of more than 
two elements and remove one of these elements. Then, we can do the same 
surgeries on |_|J points as before. This time, to recover a from a" we need 
additional step so that we increase the power of p in the bound by 2. □ 



Finally we prove Lemma A.4 by using: 



Theorem A.3 (Harer-Zagier formula [HZ86|). Let a e S 2n and define the 
genus g of a such that 

2g(a) = \a\ + |a _1 7r| — 2n + 1, 
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where ir = (1,2, . . . ,2n — 1, 2n). Then, the number of involutions without 
fixed points with genus g, denoted by e g (n), has the following recursive 
formula: 

(n + l)e 9 (n) = 2(2n - l)e 9 (n - 1) + (2n - l)(n- l)(2n - 3)e 9 _i(n - 2). 

Here, the boundary condition for the recurrence is given by eo(n) = Cat n , 
the Catalan numbers. 



Lemma A.4. We have 



e 9 (n) < 4 n ~ V 9 



for g>0. 

Proof. First, the bound is true for g = 0. This comes from the fact that 

e (l) = Cat! = 1 and (n + 1) Cat„ = 2(2n - 1) Cat„_i . 
Next, we assume that the bound is true for g — 1. Theorem [A3 implies 

e g (n) < 4e fl (n-l)+4(n-l)%_ 1 (n-2) 

< 4 [4e g (n - 2) + 4(n - 2) 2 t g _ 1 (n - 3)] + 4(n - l) 2 t g -i(n - 2) 

n-1 

< E 4 ^Vi(j-i). 

j=0 

By the induction hypothesis, 

n— 1 

e 9 (n) < ^ 1" •f V :i 
i=o 

n-1 

,3g 



i=o 



Appendix B. Measure concentration 



□ 



In this section we collect some results from asymptotic geometric analy- 
sis, which we use in this paper. 

The following theorem states that continuous functions on high-dimensional 
unit spheres have almost constant value except for sets of small measure: 

Lemma B.l (Levy's lemma HLev5110 . Let f : §> fc — > IR be a function with 
Lipschitz constant L, then there exists a constant < c < 1 and 

- Co (k-iy 



Pr{x E : \f(x) - Median[/]| > e} < exp 



L 2 
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It is well-known that singular values of random matrices of Gaussian 
entries are asymptotically concentrated [MP67J if the prime and dual spaces 
grow in sufficiently different speeds. This phenomenon is also true even 
when they are normalized to have the Hilbert-Schmidt norm one. Below, 
we quote the theorem from [ASW1 1J which was proven via Levy's lemma: 

Lemma B.2 (Aubrun, Szarek, Werner). There exist two constants < c < 
1 < C such that if we take the uniformly random unit vector \v) E C dl < 
with d 2 > C 2 d\. Then, identify \v) — V e M.di,d 2 and we have 

1 C 1 c 



Pr < { singular values ofV} (jL 



'd\ 



< exp{— cdi}. 
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